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g ! DETERMINE THE FATE OF MASSIVE STARS ? 

Cyril Georgy 1 , Sylvia Ekstrom 2 , Hideyuki Saio 3 , Georges Meynet 2 , Jose 
C ■ ' Groh 2 and Anahf Granada 2 

■ Abstract. Mass-loss rates are one of the most relevant parameters de- 

termining the evolution of massive stars. In particular, the rate at 
which the star loses mass during the red-supergiant (RSG) phase is the 
least constrained by the observations or theory. In this paper, we show 
how the mass loss during the RSG phase affects the later evolution of 
the star, as well as the final type of supernova towards which it leads. 
We also discuss some possibilities to discriminate between blue stars 
that went through a RSG phase and those which remained in the blue 
part of the Hertzsprung-Russell diagram. 
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Since more than ten years, it has become obvious that rotation, mass loss, and their 
interplay have a deep impact on the evolution of stars (see for example the review 
by Maeder & Meynet (]2012[) ). Both of them contribute to the chemical enrichment 
of the stellar surface, and the stellar winds then carry away these new chemical 
elements in the circum-stellar medium (CSM), contributing to the chemical enrich- 
ment of the Universe. The exact stellar evolution is however crucially dependent 
on several not-well-known physical processes and parametrisations, such as the 
way rotation is implemented in stellar evolution codes, the treatment of magnetic 
^ ' fields, or the mass-loss rates during the various phases of stellar evolution. 

The mass-loss rates during the red-supergiant (RSG) phase are still uncertain. 
In stellar-evolution codes, the mass-loss rates used during this phase are usually 
those determined by de Jager et al. (|1988p . Even if these rates have been confirmed 
by a recent study ( |Mauron and Josselin, 20TTj ) , some other determinations suggest 
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that in some cases, they can be much higher ( |van Loon et al., 2005 ). Observations 
also seem to show that RSGs can lose mass dramatically in some kind of bursts 
( |Humphreys et al, 2005} |Smith et al., 2009| |Moriya et al., 201 1| see also the pa- 
per by Humphreys in these proceedings). It is unclear what is the total amount 
of mass lost in the RSG phase, and how an average mass loss rate should be 
computed from the possibly discrete mass loss episodes. 

In this paper, we discuss how the mass loss during the RSG phase impacts the 
subsequent evolution of the star, and how recent stellar evolution computations fit 
(or do not fit) the observations, in the framework of single massive star evolution. 



2 Red supergiant, and then? 

2.1 How rotation and mass loss govern the duration of the RSG phase 

It has been known for long that a RSG star can evolve bluewards, if the ratio 
q = Y/° rc of the mass of its He core to the total mass becomes larger than a given 
value QGiannone, 1967[ ), typically around 0.7. In the framework of the evolution 
of a single star, rotation and mass loss affect this ratio as follows: 

• Rotation, particularly during the main sequence (MS), increases the size of 
the He core, due to the supply of fresh fuel in the core by rotational mixing 
dMeynet and Maeder, 2000p . 

• Mass loss during the RSG phase contributes obviously to the decrease of the 
total mass of the star. 

Both effects contribute to the increase of q during the RSG phase, and thus, favours 
an evolution towards higher T e g after the RSG phase. 



2.2 New models with the Geneva stellar evolution code 



Recently, we computed a new grid of stellar models at solar metallicity ( |Ekstr5m et al., 2012 \ 



with the Geneva stellar evolution code (Eggenberger et al., 2008). The free pa- 



rameters of the code (mixing length, overshooting, rotation parametrisation) were 
calibrated in order to reproduce various observations (width of the MS, red giant 
branch position, chemical enrichment of the surface of B-type stars during the 
MS). 

During the RSG phase, we use a mass-loss rate determined on the basis of 
observations by Crowther (|2001j) . which is similar to the rate of de Jager et al. 
(1988). However, for the most massive stars of this grid (Mj n j > 20 Mq), we 
observed that some layers in the outermost part of our models became significantly 
more luminous than the Eddington luminosity of the star, due to the strong opacity 
in the iron peak. In that case, this very unstable situation could favour mass 
loss. We thus decided to increase the standard mass-loss rate by a factor of 3 in 
that case. Interestingly enough, it led to mass-loss rates comparable to the ones 
determined by van Loon et al. ()2005j) . This kind of increased mass-loss rates 
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Fig. 1. Hertzsprung-Russell diagram tracks for rotating models (Z@) between 5 and 
40 Mq (Ekstrom et al., 2012). The red dots correspond to the position of Galactic RSGs 
(Levesque et al., 2005) 



during the RSG phase was also used in previous works ( |Vanbeveren et al., 1998a| 
|Vanbeveren et al., 1998b[ ), and more recently in Chieffi & Limongi (j2013p . that led 
to similar conclusions. 



2.2.1 Consequences on RSGs 

The inclusion of an increased mass-loss rate for the more massive RSGs at solar 
metallicity has as the direct consequence to help these models to evolve blue- 
wards at some points during their RSG phase. On Fig. [T] are shown the tracks 
in the Hertzsprung-Russell diagram (HRD) for our rotating models between 5 
and 40 Mq. The red dots correspond to the positions of observed Galactic RSGs 
(Levesqu e et al., 2005[ ). First, we see that our models reproduce well the observa- 
tions, particularly the upper luminosity of the RSGs. 

The second interesting point is that rotating models of 20 and 25 Mq evolve 
away from the RSG branch at the end of their life. This is a direct consequence 
of the strong mass loss they encounter during the RSG stage. 
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2.2.2 Blue-supergiants and Wolf-Rayet stars, and consequences on SN progeni- 
tors 

From our grid of rotating models, the following evolutionary scenarios are possible 



(Georgy et al., 2012): 



• MS -» RSG for star with initial mass M < 16.8 Mq 

• MS -> RSG -> nitrogen rich (WN) Wolf-Rayet star (WR) for star with initial 
mass 16.8 M < M < 25.0 Mq 

• MS ^ Yellow supergiant (YSG) -)• WN ->■ carbon rich (WC) WR star for 
star with initial mass 25.8 Mq <M< 60.0 M q 



MS WN ->■ WC for star with initial mass M > 60.0 M, 







Interestingly, the models between 20 and 25 Mq end their life as WN stars 
(or even as luminous blue variable (LBV) stars, see Groh et al. (|2013|) V with 
log(T c ff) ~ 4.1 — 4.3. This is in good agreement with the identified progenitor 



of SN 2008ax ( Crockett et al., 2008 j). The estimated m aximal mass we obtain for 
a type IIP progenitor of 16.8 Mq ( |Georgy et al., 2012 1 is also in good agreement 



with the work of Smartt et al. (2O0S]), who determine a maximal initial mass of 
16.5 ± 1.5 Mq for the progenitors of this kind of Supernovae (SNe) (Note however 
that the so-called "red-supergiant problem" could have other explanations, see 
Walmswell & Eldridge ([2UT2"[n . 

2.3 Strong mass-loss rates for less massive RSGs? 

During the past 20 years, 3 progenitor stars have been identified for some (rare) 
type IIL or lib S Ne: SN 1993 J (|Maund et al., 2004| ), SN 2009kr QElias-Rosa et al., 2010D , 
and SN 2011dh ( |Maund et al., 20fT| ), as well as a peculiar progenitor for the type 
IIP SN 2008cn ( |Elias-Rosa et al, 2009] ) . All these progenitors are YSGs, with 
3.63 < log(T ff) < 3.78, and a luminosity of log(i/L Q ) ~ 5. The luminosity indi- 
cates a star with an initial mass of ~ 15 Mq. However, single star models show 
that with such a mass, the star ends its life as a type IIP SN. Binarity could play 
an important role in the evolution of such stars, as the progenitor of SN 1993J 
is most probably a binary star QPodsiadlowski et al., 1993] ), and this possibility 
is also not excluded for SNe 2008cn and 2009kr. However, the observations for 



SN 201 ldh do not indicate any hint of binarity (Maund et al., 2011). A recent 
work by Benvenuto et al. (12013[ ) shows that the secondary star could be hidden, 
thus the absence of a companion cannot be taken as a proof that the star did not 
interact with a companion in the past. 

We can however wonder if a single star model, with an increased mass-loss rate 
during the RSG stage, is also able to reproduce the position and characteristics 



of these progenitors (see ( Bersten et al., 2012] ) ). In Georgy (I2012p . we have shown 



that it is indeed the case. Fig [5] shows the tracks in the HRD for a 15 Mq model 
at solar metallicity, with a standard mass-loss rate (red solid track), or with a 
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Fig. 2. HRD tracks for 15 Mq rotating models (Zq, Georgy, 2012) with various mass-loss 
rates during the RSG phase: standard one (red solid curve), 3 times (green dot-dashed 
curve) , 5 times (purple short-dashed curve) and 10 times standard one (blue long-dashed 
curve). The observational data are from Maund et al. (2004, 2011), and Elias-Rosa et 
al. (2009, 2010). 



mass- loss rate increased by a factor of 3 (green dot-dashed track) , 5 (purple short- 
dashed track) or 10 (blue long-dashed track), during the RSG phase. We see that 
the effect of increasing the mass-loss rate is to shift the final position of the model 
towards higher and higher T e g. The final properties of the progenitor star (core 
mass, ejected H mass) are also compatible with the values deduced from the SN 
light curve QBersten et al., 20120 . 

Note also here that increasing the mass-loss rate during the RSG phase shortens 
this phase, and thus the duration of the episode with strong mass loss. These effects 
compensate somehow, and the models with 3 to 10 times the standard M finish 
with roughly the same final mass. Our computations suggest that this effect forbids 
a later evolution to a WC star for these lower-mass RSGs (see Meynet et al. in these 
proceedings). This last point is however sensitive to the detailed mass- loss his- 
tory, as some other authors QVanbeveren et al., 1998a| |Vanbeveren et al.~9 98b) 
obtained such WR stars for models with initial mass as low as 15 Mq. 
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Table 1. Main abundances before and after the RSG phase of a non-rotating and a 
rotating 25 Mq model. By "before" and "after" , we mean the first and last time that 
log (T ff) < 3.66 (roughly corresponding to the temperature at which the surface convec- 
tive zone begins to sink). 



element 


non rotating 
before RSG after RSG 


rotating 
before RSG after RSG 




0.720 


0.450 


0.638 


0.370 


4 He 


0.266 


0.536 


0.348 


0.616 


12 C 


2.28- 10" 3 


5.52- 10" 5 


1.07- 10~ 3 


1.97- 10~ 4 


14 N 


6.59 • 10" 4 


7.32 ■ 10" 3 


3.54 • 10~ 3 


6.51 ■ 10^3 


16Q 


5.72 • 10~ 3 


1.10- 10~ 3 


3.94 • lO-' 3 


1.81 ■ 10^3 


12 C/ 13 C 


81.43 


3.58 


9.47 


5.42 



2.4 How to distinguish a blue star coming from the red and an "ever-blue" 
star? 

Is there a way to determine whether a blue star evolved from a previous RSG 
phase? For example, is it possible to determine observational properties that are 
unique to a star coming from the red part of the HRD compared to another one, 
which has still not crossed it? 

A first point to notice is that the surface chemical composition would be dif- 
ferent. A star which has been previously in a RSG phase has encountered the 
first dredge up due to the large convective zone developing during this phase. The 
effect is to bring to the stellar surface products of H-burning through CNO cycle, 
mainly, a strong enrichment of He and N, and a depletion in H, C and O. Ta- 
bled] shows the abundances of these elements, as well as the ratio 12 C/ 13 C for the 
non-rotating and rotating 25 M Q models of Ekstrom et al. (|2012p . 

Another indication could come from the pulsational properties of the star. On 
Fig. |31 we indicate by red dots models which have at least one radial pulsational 
mode, for four different initial masses. An interesting thing is that during the first 
crossing of the HRD, no modes are excited, and the star should not vary. However, 
several a Cyg variables are found in the corresponding T c g region. On the contrary, 
if the star crosses again the HRD towards the blue side, the higher M/L ratio (due 
to the strong mass loss during the RSG phase) allows some pulsational modes to 
be excited in the T e g range between ~ 3.7 and ~ 4.3. It could thus be possible 
to discriminate between these two types of stars on the basis of their pulsational 
properties. This will be studied in detail in a future work (Saio et al., submitted). 

Finally, due to the strong breaking encountered during the RSG phase, stars 
coming back from the red are probably slow rotators, contrarily to the stars which 
have still not crossed the HRD, which can have any rotational velocities. 
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Fig. 3. HRD tracks for various models (non-rotating ones in dotted lines, rotating ones 
in solid lines). The 14 Mq model was computed with a mass-loss rate increased by a 
factor of 5 during the RSG phase. The red dots indicate models for which at least one 
radial pulsation mode is excited (Saio et al., submitted). The blue points are observed 
q Cyg variables in the Galaxy (Kudritzki et al., 1999; Przybilla et al., 2010; Firnstein 
and Przybilla, 2012) and in NGC 300 (Bresolin et al., 2004; Kudritzki et al., 2008). Blue 
line indicates /3 Cephei instability region (log(T c ft) ~ 4.3 — 4.5) and the blue-edge of the 
Cepheid instability strip (log(T'eff) ~ 3.8). 



3 Conclusions 

The mass lost during the RSG phase is a key parameter determining the later evo- 
lution of massive stars going through such a phase. We have seen that increasing 
the standard mass-loss rates during the RSG phase for the most massive RSGs 
allows for a good fit of the maximal luminosity of Galactic RSGs, and also of the 
maximal initial mass leading to a type IIP SN. At lower mass, it also provides an 
alternative scenario to the binary channel to explain the peculiar position in the 
HRD of the progenitors of type IIL or lib SNe. 

Finally, we also propose some observational tests which could make it possible 
to distinguish the blue stars that evolve from a previous RSG phase and other 
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ones, that did not have a previous RSG stage. 
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